RIGOROUS DERIVATION OF THE HYDRO DYNAMICAL 
EQUATIONS FOR ROTATING SUPERFLUIDS 



HAILIANG LIU AND CHRISTOF SPARBER 

Abstract. Using a modified WKB approach, we present a rigorous semi- 
classical analysis for solutions of nonlinear Schrodinger equations with rota- 
tional forcing. This yields a rigorous justification for the hydrodynamical sys- 
tem of rotating supcrfluids. In particular it is shown that global-in-time semi- 
classical convergence holds whenever the limiting hydrodynamical system has 
global smooth solutions and we also discuss the semi-classical dynamics of sev- 
eral physical quantities describing rotating supcrfluids. 



1. Physical Motivation 

Bose-Einstein condensates (BECs) play an outstanding role in present-day physics, 
cf. [Hj for a general introduction. Understanding and controlling their behavior 
is of great fundamental importance and essential for novel applications. A partic- 
ular focus of interest is on dynamical phenomena related to superfluidity and the 
creation of quantized vortices, see, e.g., [T| 112 ] [15 ] 116]. To this end the typical ex- 
perimental set-up is based on a trapping potential subject to a rotational forcing. 
The mathematical description of this system is then usually given by the celebrated 
Gross-Pitaevskii equation (GPE), a mean field model for the macroscopic wave func- 
tion of the condensate (see [17] for a rigorous derivation). In the rotating reference 
frame we are consequently led to 

h 2 

(1.1) ihdtip = - — + V (x)i> + k#| 2 V> + in x x • V X V, 

with a nonlinear coupling constant k = NAnha/m, where N is the number of par- 
ticles forming the condensate and m, a e R respectively denote their corresponding 
mass and scattering length. In the context of BECs the potential V is usually as- 
sumed to be a harmonic oscillator confinement, i.e. Vq(x) = ^mujQ\x\ 2 , for some 
ujq G M. d . In p.ip we also write x L = (x2, — xi,0) T in d = 3 spatial dimensions 
and analogously x 1 - — (x2, — £i) T for d = 2. The latter case is motivated by recent 
experiments for BECs which are strongly confined in one or two directions. Thus 
x ■ Vz can be interpreted as the negative ^-component of the quantum mechanical 
angular momentum operator L = x X (— iV^) and f^o is the corresponding angular 
velocity. 

In the so-called Thomas-Fermi limit, which typically applies for large systems, the 
dynamics is presumably well described by the hydrodynamical equations for rotating 
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superfluids [HHIM], i.e. 

d t p + V x • (p(v - fioar 1 -)) =0, 

9tv + f4r- - • v + y + P ) = 0, 



(1.2) 



where p := |^| 2 denotes the particle density and v the corresponding supcrfluid 
velocity defined by 

K ImfiJVxV) 
v := — . 

m |i/)| 2 

The passage from (|1.1|) to (|1.2|) is usually explained by using the classical Madelung 
transformation of the wave function, where one writes 



(1.3) #,!) = #, x) exp (i<I>(i, , 

and consequently identifies v := V 2 -$. Formally plugging the ansatz (| 1 . 3j) into (|1.1|) . 
separating real and imaginary parts, and discarding terms oc H 2 yields (|1.2p . This 
asymptotic regime is particularly interesting for numerical simulations, cf. [HG3, and 
we also refer to |18[ lllj for an extensive review of such dispersive limits. However, 
the representation (|1.3p makes no sense at vacuum points, where p = 0. Indeed the 
system (|1.2I) degenerates at such points and is therefore only weakly hyperbolic (see, 
e.g., HH]). 

In the present work we shall rigorously prove that (|1.7p approximates (|1.1[) in a 
certain sense and moreover draw some further conclusions from that. To this end 
we rescale equation (11.1[) , as described in [5l [20] , into its dimensionless form. This 
yields 

,2 

(1.4) \edti> E = -— + V(x)^ e + eS 5/2 \ip e \ 2 ^ +knx- L • V.f, 

with e = H/(mL 2 ) and S = A-naN/a®, where L denotes the characteristic length 
of the condensate and ao = y/H/(muio) is the ground state length of the harmonic 
oscillator potential V(x). From now on and V(x) denote (rescaled) dimension- 
less quantities. The particular choice |e(5 5 / 2 | = 0(1), which yields ^ « L and 
thus e = (ao/47r|a|A f ) 1 / 5 <C 1, corresponds to the Thomas-Fermi regime of strong- 
interactions, cf. [5] for more details. Note however that our scaling is different from 
the one used in [TJ [13) . 

We consequently study the following semi- classically scaled nonlinear Schrodinger 
equation (NLS) with rotational forcing 

(15) Ued t r= -^Ar+v( X )r+f(\r\ 2 w+^Qx ± - v x r 

where ( e 1, i £ R d , for d = 2, 3, and f2 > 0, some given constant (independent of 
Kl). The energy functional associated to (|1.5p reads 

(1.6) E(r)= [ ^\V x ^\ 2 + {V(x)+F(\r\ 2 ))m 2 + ^nWx ± -V x r dx, 

JR d * 

where F is the primitive of /. From now on we impose the following assumptions. 
Assumption 1. It holds: 

• The nonlinearity satisfies f G C°°(R) such that f > 0. 

• The potential V is of the form V(x) = ||o; • x\ 2 , with ui G M. d . 
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The supcrfluid equations corresponding to the NLS p.5[) arc thus given by 
f d t p + V x - (pCv-nar 1 )) =0, 

( L7 ) S / Ivl 2 \ 

I d t v + V x (LL-nx^-v + V + fip)) =0. 

We expect the system (|1.7[) to be a valid description of solutions to p.5[) in the 
limit e — ► 0. In order to prove this rigorously, we will heavily rely on semi-classical 
expansion techniques first developed by Grenier |19j and later extended by Carles 
in [10] . The latter work in particular treats the case of a harmonic confinement 
and nonlinearities of different strength (see also [3] for an extension to higher order 
nonlinearities and [14] for results on the derivative NLS). We can thus focus on the 
role played by the rotational forcing. Indeed the motivation for our work is threefold: 
First, we aim to generalize the results of [10] to the case where a rotational term is 
included and consequently give a rigorous justification of (jTTTJ) . Second we strengthen 
these results a bit in the sense of Corollary 13.11 below. Finally we aim to describe 
the dynamical features of rotational BECs from the semi-classical point of view, as 
given in Theorem 14. II 



2. The modified WKB-approach 



As noticed in [TO] [19] , the classical Madelung transformation (|1.3p is not well suited 
to rigorously derive the semi-classical asymptotics of ip e . Rather one is led to consider 
a modified version of it. To this end one writes the exact solution to (|1.5p in the 
form 

(2.1) ip e (t,x) = a e {t,x)e^^ x)/e , 

where from now on the "amplitude" a e is allowed to be complex-valued. Moreover a e 
as well as the (real-valued) phase $ e are assumed to admit an asymptotic expansion 
of the form 

(2.2) a e ~ a + eai + e 2 a 2 H , $ e ~ $ + e$i + e 2 $ 2 H . 

Since a e is complex valued, the phase $ £ can be seen as an additional degree of 
freedom introduced as a multiple scales representation for ip e . In any case, the 
ansatz (|2.ip should not be confused with (|1.3[) and in particular it has nothing to 
do with a rewriting of ip e into polar coordinates. The main gain of this modified 
WKB-approach is that it yields a separation of scales within the appearing fast, i.e. 
e-oscillatory, phases and slowly varying phases, which eventually can be included in 
the (complex-valued) amplitudes. 

An analysis based on WKB-type methods necessary requires admissible initial data. 
To this end we introduce the following definition: A function / <G C°°(IR rf ) is said to 
be sub- quadratic, if 

(2.3) VaeN d , H >2: sup \d«f\ < oo. 
We consequently impose: 

Assumption 2. The initial data is of the form 

(2.4) ^ n (x) = at n {x)e^^. 

Here af n is complex-valued and admits an asymptotic expansion in e whereas < f>; n (x) 
is e-independent, real-valued, and sub-quadratic. 
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It is important to note that Assumption[5]imposes a particular e-oscillatory structure 
on the initial data but does not yield any problems at vacuum points since ct? n is 
allowed to be complex-valued. 

Remark 2.1. In particular we are free to choose, say in d = 2 spatial dimensions, 
an initial data of the form: $i n (x) = 0, af n (- T ) = x( r ) eimS : where we have used 
polar coordinates. Here m £ No is the so-called winding number. For m ^ and 
X appropriately chosen (see, e.g., |H[H1[S]) this allows for so-called vortex states as 
initial data. We also note that we could allow for more general initial phase function 
3>; n which admit an asymptotic expansion in powers of e (analogously to af n ). For 
simplicity we do nit treat this case though. 



Upon substituting (|2.ip into the NLS (|1.5p . we have the freedom to split the appear- 
ing terms into 




\\v x &\ 2 + v(x) - ttx 1 - ■ v x ^ + /(K| 2 ) = o, 

V x $ e • V x a e + ia e A$ e - tlx 1 - ■ V x a e = ^Aa e . 



(2.5) 



This system is equivalent to the nonlinear Schrodinger equation (|1.5|) . In particular 
it does no longer represent a splitting into real and imaginary parts (since a e is 
complex- valued), in contrast Madelung's original approach. Moreover the system 
(|2.5[) is seen to be perturbed by a term which is "only" 0(e) (since the term oc Aa 
now appears in the equation for the amplitude instead of the one for the the phase). 

As e — > 0, the corresponding limiting WKB system is then (formally) found to be 



(2.6) 




-Six*- ■V x <£ + f(\a\ 2 ) = Q, 
aA$ - fix 1 - • V T a = 0. 



From here, setting p = \a\ 2 and v ~ Vj;$, one (again formally) obtains a hydrody- 
namical equations of superfluids (|1.7p . Indeed we shall prove below that (|2.6p and 
(|1.7p arc in a certain sense equivalent. 



Remark 2.2. From the point of view of geometrical optics, the above given limit 
corresponds to the supercritical case, cf. [TQ1 [SI LH] for more details. Note however 
that a rotational forcing is neither considered in [TU1 [3] , nor in [TH] . Also note that 
in the case of a linear Schrodinger equation the corresponding system (|2.6p would 
be decoupled since the first equation would be replaced by the classical rotational 
Hamilton- J acobi equation (HJ) 

(2.7) d t S + ^\X7 X S\ 2 + V{x) - tlx 1 - ■ V X S = 0. 

Equations of the form (|2.7p have been extensively studied in 20J, where several 
qualitative properties for the corresponding solutions are established. 

As a first step in our analysis we shall show that the usual hydrodynamical system 
(|1.7p and the WKB system (|2.6p admit smooth solutions on the same time-intervall. 

Lemma 2.1 (Equivalence). Let T be the maximal time of existence for a smooth 
solution (/O, v), with p > 0, of the hydrodynamical system (jTTTJ) and let T* be the 
maximal existence-time of a smooth solution (a, $) of (|2.6p . Then we have T = T* . 
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Proof. We define (a,(3,v) := (Rea,Ima, V x $) and rewrite (|2.6[) in the following 
form 



(2- 



' d t a + v ■ V x a + ^V x -v - fix- 1 - • \7 x a = 0, 
filt/9 + v • V x f3 + | V x • v - fix x • V x /3 = 0, 
9 t v + V, f ^ - Qx 1 - ■ v + V + f(a 2 + /3 2 ) ) 



Assume that (a, /?, v) is a smooth solution of the modified WKB system (|2.8[) for 
t £ [0,T*). Then r] = a 2 + (3 2 and v = are smooth solutions of (JTTTJ). By 

uniqueness it follows that (77, v) = (p, v) for t < T* and hence T* < T. 

Conversely, assume that (p, v) is the smooth solution of the hydrodynamic system 
(|1.7[) for t e [0, T], subject to initial data {p- m > 0, v in ) such that p; n = a? n + (3f n . 
By assumption the velocity v is smooth, one thus obtains smooth a and b from the 
first two transport equations in (|2.8p . subject to initial conditions (ai n ,/3i n ) T . A 
combination of the two equations for a and (3 gives 

d t (a 2 + (3 2 ) + v V^a 2 + /3 2 ) = -(a 2 + /3 2 )V X • v + fix 1 - ■ V x {a 2 + {?). 

Subtracting this from the equation for p in (|1.7p . we find that p := p — (a 2 + (3 2 ) 
solves a transport equation 

d t p + • (v/5) = Oa;- 1 • V x p , 

with /5| t = 0, and hence p(t, x) = 0, i.e. p = a 2 + (3 2 . This shows that (p, v) is also 
the smooth solution of the hydrodynamic system (|1.7|) for i € [0, T*], hence T < T* . 

In summary this yields T = T* for the solutions of (|2.6|) and (|2.8|) . To get back to 
$ itself, we first note that in the equation 

+ ^IV^I 2 + V(x) - tlx 1 - • + f(\a\ 2 ) = 

all terms are uniquely determined by (|2.8j) except for <9t$. Imposing 3?L_ = 3>i n , 
such that Vi n = V^^in, and setting 

$(t,a:) = *to{x)+ J Q|v(r,:r)| 2 + \a(r,x)\ 2 + V(x) - Qx x ■ v{r,x)j dr, 

we infer that d t {V x $ — v) = V x 9 t <I> - v = 0, hence v = V K <i>. This then fully 
determines $ on [0,T]. □ 

Having established the equivalence between the limiting systems (|1.7p and (12. 6p . 
we can now focus on deriving rigorously (|2.6p from (|2.5[) (or, equivalently, the NLS 

(HD). 

3. The classical limit 

From now on we shall mainly consider the WKB system (|2.6[) in the form (|2.8p which 
allows for a treatment in the sense of hyperbolic systems. In the corresponding 
analysis, the rotational HJ equation (|2.7|) becomes important. As a preparatory 
step we shall therefore study the Cauchy problem 

f 3t5 + ^IV^I 2 + 7(x) - Slx x ■ V X S = 0, 

I ^lt=0 = *-'in(2 ; ), 

The following result will be used throughout this work. 
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Lemma 3.1. Let Assumption^ hold. If S in (x) G C°°(R d ) is sub- quadratic, then 
there exists a r > such that \2. 7| ) admits a unique smooth solution for t G [0,r). 
Moreover, the phase S(t, x) remains sub-quadratic in x, for all t G [0, r). 

As pointed out in [lOj , the sub-quadratic assumption for the initial phase is sharp for 
solving (|2.7[) globally in space even in the presence of no rotational force. Concerning 
the existence of global- in-time smooth solutions we refer to [20] . 

Proof. The Hamiltonian corresponding to (|3.1j) is 

H = ^\p\ 2 + V(x) - Oar 1 - p. 
and thus the corresponding Hamiltonian flow is governed by 

jx = p-flx 1 -, x\ t=0 = x , 

\ p = -uix - flp x , p\ t=Q = V a 5i n (x ). 

Standard ODE theory implies that there exists a unique solution for t G [0, r], in 
which the map x = x(t, xq) is well defined, satisfying 

d 

dct(r) > 0, r := — x(t,x ). 

dx 

Existence of smooth solutions is therefore guaranteed and we now have to prove that 
S remains sub-quadratic for t G [0,r]. From [20] we infer that the phase gradient 
u = \7 X S satisfies 

dtu + (u — fix 1 ') ■ V x u — —ljx — flu 1 ' . 

Further, along the particle path induced via x = u — fix 1 - the Hessian of the phase 
function £ := D^S(t, x) solves a matrix ODE 

(3.2) AS + S 2 = -ujl d , 

where we shortly denote 

D t :=d t + (U-Qx x )-V x . 

This shows that S remains uniformly bounded in terms of xq G M. d . Finally we show 
that the existence time r does not shrink as x varies over M. d . Differentiation of 
x = u — fix 1 - in terms of xq yields 

D t r = (s - nj)r, j = 

for d — 3. (The expression for J in the case d — 2 is obvious.) This gives 

dct(r) = cxp (J^ Tr(E - flJ) dr^j = exp (J^ Tr S dr j , 

where Tr denotes the standard trace map. The fact that det(r) does not depend on 
Xq explicitly implies that S is uniformly bounded in xq G K d . □ 




In order to prove the local- in-time existence for the system (|2.5[) . we follow the 
strategy in [10] and decompose the phase <3? £ into 

(3.3) $ e = ip e + S, 
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where S is the smooth, sub-quadratic phase function guaranteed by Lemma |3. II In 
terms of ip e and a e the system (|2.5|) becomes 



(3.4) 



d t <p e + v x s ■ V x <p e + i|V^ e | 2 - flx x • V x <p e + /(|a £ | 2 ) = 0, 

d t a c + V X (S + ip e ) ■ V x a e + %- A(S + ip e ) - fix 1 - ■ V x a e = ^ A x a e . 



Note that this set of equations is still equivalent to the nonlinear Schrodinger equa- 
tion (|1.5|) . The reason for decomposing $ e via (|3.3[) is rather technical and due to 
the inclusion of the potential V and the rotational term. 

For notational convenience we further introduce 

W := (Rea £ , Imo £ , d Xl tp e , . . . , d Xd ^) T , 

where a £ and y £ satisfy (|3 .4[) (an analogous notation is used for the corresponding 
initial data). Moreover, we shall frequently use the notation 

(3.5) mU € {t])--=\\U e {t)\\ s + \\\AU'{t)\\ s -v 

where || • || s is the usual i/ s (R d )-norm. 

Proposition 3.1 (Local existence). Denote by r > the existence time of smooth 
solution S(t,x) to (|3.ip . and let the Assumption]]] hold. Consider the Schrodinger 
equation (|1.5p subject to initial data, which satisfy Assumption [D such that U* n G 
H s (R d ) and \x\Uf n G H s - 1 (R d ), for s > 2+d/2. Then there exists a time T £ G (0, t), 
and a unique solution to (|1.5[) of the following form 

ip e {t,x) = a e (t,x)e i ' s '' {t ' x)/e , for 0<t<T e . 

Moreover, it holds 

U e G L°°((0,T e ],H s (R d )), \x\U e G L°°{{%T^,H a - l (pL d )). 

We thus know that locally-in-time the oscillatory structure of the modified WKB 
representation for solutions to Q1.5P persists, as long as the classical rotational HJ 
equation has smooth solutions. 

Proof. Introducing the velocities v e := W x ip e and w := S7 X S — fix 1 - in (|3.4[) . we have 

C d t v e + {v e + w) ■ \7 xV e + {V x w)v e + V x /(|a £ | 2 ) = 0, 

(3-6) < a e ie 

{ d t a e + (v" + w) • V x a £ + — • (w + v e ) = -A x a e . 

From Lemma 13.11 we know that w = V^S* — fix 1 - is indeed sub- linear. We further 
separate a £ into its real and imaginary part, i.e. a £ = a £ + i/3 £ , to obtain the following 
hyperbolic system 

d 

(3.7) d t W + J2 iMU") + B 3 (w)) d X] U £ + M(V x w)U e = | LU e , 

3=1 

where U e = (a e , [3 e , v\ , • • • , v d ) T . The coefficients matrices are 

d 

J=l 
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with /' = /' (|a £ | 2 + |6 e | 2 ) and 

/ -w 
M(V x w) = \V X ■ w | , L = 

\ V x v 

Observe that (|3.7[) can be symmetrized by 




^Lj 





1 

4/' 

which explains the necessity for our assumption /' > 0. (In 3\ it is shown how to 
overcome this difficulty in case of higher order nonlinearities.) We are thus able to 
proceed with energy estimates in the Sobolev space H s (R d ), which follow from the 
classical theory for hyperbolic systems: 

Lemma 3.2. Assume that w is sub-linear. For s > 2 + d/2, \a\ < s and \(3\ < s — 1, 
there exists a locally bounded map €{■), satisfying C > and C(0) > 1, such that 

(3.8) 1 (QffZU',d2U*) <C(N[U<}) \\u*\\ 2 s , 

(3.9) ^ t (Qde{ Xj U%dl{ Xj W)) <C{N[U*]) (irt + UMEHI^) . 
where (■, •) denotes the usual scalar product on L 2 (M. d ). 

Equipped with these estimates, we are able to conclude the local-in-time existence 
result. Set 

d 

\a\<s j=l \0\<s-l 

then the estimates presented above yield 

£[U e (t)} <£[U e (0)} + C f C(M[U e (T)})£[U e (T)}dT. 
Jo 

Invoking a Gronwall-type inequality we consequently arrive at 

£[U e }{t) < £[U C {0)} exp (c C{Af[U e ]{T))dT 

which we expect to have a relaxed bound 2£(0) for a finite time where J\f[U e ] < C\ 
for C\ > Co = Af[U* n ]. Thus an existence time exists and satisfies 

In 2 

t < T f 



CC(d) 

It is obvious that the smaller the initial data (measured by Ci), the larger the 
time-interval of existence. A local-in-time existence is thus established. □ 

Remark 3.1. The local- in-time existence for solutions of the limiting system (|2.8[) 
can be proved analogously. 

For e G (0,1], assume that the initial data satisfies N[U* n ] < Co < C. Thus, for 
e G (0,1] fixed, the local existence shows that for any number C\ G (Co,C), there 
exists a T e > so that (|2.6| has a unique classical solution satisfying J\f[U e ] < C\ 
for t G [0, T e ]. Define 

T e := sup{0 < T e < t : Af[U e ]{t) < C u Vt G [0,T e ]}. 

Namely, [0, T e ) is the maximal time- interval of existence and depends on C±. It will 
be necessary to show that lim c ^o2 n£ > 0, which we shall prove in Corollary 13. II 
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Theorem 3.1 (Convergence rates). Under the same assumptions as in Proposition 
lff.il suppose that there exist ai n , $i n £ H s (M. d ), for s > 2 + d/2, such that 

\\at n -a in \\ s = 0(e). 

Let U := (Re a, Lma, d Xl ip, . . . , ^ip) 1 be the smooth solution to ()2.6p corresponding 
to the initial data (4>i n ,ai n ). i/ 

17 e L°°({0,T*],H s (R d )), \x\U e L oo ((0,T*],H s - 1 (R d )) 

wi//i T* > finite, then there exists eo and C» > 0, smc/i £/ia£ /or e < eo 

II a e W - a(t) || s < C„e, || - S(t) || s < Cet, 

/or a/a 6 [0,min{T*,T e }). 

As a consequence, we infer that the pair (p,v) = (|a| 2 ,V x 4>), solves the hydrody- 
namical system (|1.7[) . Note that $ = <p + S, where 5 is determined by the HJ 
equation (|2.7p . 



Proof. As C/ e solves the forced hyperbolic system p.7[) , then the corresponding lim- 
iting function U is governed by 

d 

(3.10) d t U + ^(A,(/7) + Bj(w)) d Xj U + M(V x w)U = 0. 

i=i 

Denote l¥ e := U € - U, we thus have LU € = LW £ + LU and W £ solves 
d 

(3.11) d t W £ + J2( A A u<i ) + Bj{w))d Xj W e + M(V x w)W e = | LW £ + Ft 

3=1 



where 

d 

R £ := ^ (^-(t0 - ^ ([/<)) 9 a ,Ci + | L£7 
i=i 

Using Lemma [3~2l we obtain, for s > 2 + d/2, 

(3.12) A(Qc>«W e ,^Vn <C(iV[^])||^ e ||2 + (g^ J R e ,^^), |a|<s, 
at 

as well as 

(3 13) ^ (*^),df ( Xj W*)) < C(N[W £ ]) (\\W £ \\ 2 S + || \x\W £ \\U) 

+ (Qd^(x j R%d^x j W e )) , 

for |/3 1 < s — 1. Now we estimate terms involving i? e in (|3.12p and (|3 . 13() . This is 
done only for t G [0, min{T*, T £ }), in which both U e and U are regular with uniform 
bounds for Af[U £ ] and AT[{7]. 

First we have 

{QdgB?,d2W*) < C\\d%W% \\d%R% 

and 

||S^|| < C||^C/|| s ||^(C/ e )-^([i)|| w +Ce||^C/|| . 

Note that 

A J {U £ )-A 3 {U) = {W £ ) 3+2 \ d+2 + A, 
where the only non-zero entries of A are A\ t j+2 = \{W £ )i, A 2 .j+ 2 = 2 -(W £ ) 2 and 
A m = 2f'(\a £ \ 2 )a £ - 2f'(\a\ 2 )a, A j+2 , 2 = 2/'(|a e | 2 )6 e - 2/'(|a| 2 )6. 
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Using these relations and the boundedness of ||£/|| s+ i we conclude that 
\\Ai(U c ) - MU)\\ la] < C(M[W% \\U\\ S+1 ) \\W*\\ S . 

Therefore 

^(Q^S9^)<C(Ar[^]J|^|| s+ 0||^lf + Ce||C/|U|^|| s . 

|a|<s 

Further calculations give 

]T (Qdgix j IF),dg(x j W e )) <C(N[W% \\U\\ S ) (\\W% + \\\x\W'fc_^ 

\0\<S-1 

+ Ce\\\x\W%_ 1 . 

Substituting all these estimates into (|3. 12|) and (|3 . 13[) and integrating over M. d we 
obtain 

£[W e (t)}<£[W e (0)} + C f C{N[W%\\U\\ s+1 )£[W e {T)]dT + Ct f M[W e (r)] dr, 

Jo Jo 

where we have used the fact that £ [W e ] and A/" 2 [W^ e ] are equivalent in the sense that 
there exists a constant C > such that C^Af 2 ^} < £[W £ ] < CAf 2 [W e }. Note 
that £[W £ (0)} = 0(e 2 ), we thus have for t £ [0, T* A T] 

£[W e (t)] <C(T*)e 2 +c f C(N[W%\\U\\ s+1 )£[W'{t]&t. 
Jo 

We apply Gronwall's inequality to obtain 

£[W e (t)]<Ce 2 cxp (cj C{N[W%\\U\\ s+1 )£[W e {T)}AT^j =: T(t). 
Thus, we infer 

±F(t) <CW{W*](t),\\U\\ s+1 (t))F(t) 2 , T(0) = Ce 2 . 

This differential inequality with 0(e 2 ) initial data ensures that for fixed T*, there 
exists e and C* such that for t G [0, T*} and e < eo we have 

£[W e (t)] < T(t) < (C*e) 2 . 

This leads to 

|| a\t) a(t) I = 0(e), \\V x ^(t) V x *(t)|| a = 0(e). 
The second estimate combined with the equations for $ e and <f> yields 

wnt)-m\\ s =o(e)t 

and the proof is complete. □ 

Theorem 13.11 yields an approximation of "0 e for small times only |10j . since 
|| r(t) - a(t)j*®' e \\ L2 = || fl £ (i)e*( f '/ £ - a(t)j*W*\\ L2 

< || a*(t) - a(t) \\ L2 + || e^V* - J*W<\\ La . \\ a(t) \\ L2 

< 0(e) + 0(f). 

In other words, to accurately approximate the wave function tp e itself one has to 
take into account higher order corrections. Indeed, it has been shown in |10j that 

(3.14) || r(t) - «(t)e i4l(t) e W(i)/s L-ao^L'm) < 0(e), 

where $i is the first corrector appearing in the asymptotic expansion (|2.2p . The 
slowly varying, phase <&i(t,x), for t G [0, T*], is obtained from a linear hyperbolic 
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system with source terms (and thus will always be generated during the course of 
time). Consequently one might consider a := ac 1 * 1 as a new (complex- valued) WKB 
amplitude. Since in the present work we are mainly interested in deriving (|1.7[) . we 
shall not go into further details and rather refer to (TOl [3] . Note however, that in 
the case where one aims to accurately describe the semi-classical dynamics of vortex 
states (see Remark 12. ip it is crucial to take into account this additional slowly 
varying phase <£>i. Even though the fluid velocity in the classical limit is given by 
as usually, the total phase of the wave function in this asymptotic regime is 
4>-|-e<I>i . This 0(e)-correction of the phase is usually ignored in the physics literature. 
However, by doing so one can no longer justify a semi-classical approximation in the 
sense of (|5TI)) . 

With the above given result in hand, we are now able to show global-in-time con- 
vergence of the semi-classical limit. 

Corollary 3.1 (Global convergence). Under the same assumptions as before and 
for any C\ satisfying 

(3.15) N[W ] < C < Ci, N[U e (t)] < Ci < C, for t G [0, mm{T*, T £ }), 
it holds T e (Ci) > T* for e > sufficiently small. 

In case the hydrodynamic system ()1.7|) is proved to admit global solution, i.e. 
T* = oo, we have consequently established convergence of solutions of (|2.5|) towards 
solutions of (|1.7|) globally in time. 

Proof. Assume the contrary of what we aim to prove, i.e. assume that there is a C\ 
satisfying (|3 . 1 5|) and a sequence e„ — > as n — » oo and T e ™(Ci) < T*. Then there 
exists C satisfying 

N[U{t)] < C < d. 

From Theorem 13.11 it follows that 

M[U f -{t) - U{t)] < CU. 

Thus, there is a n G N such that Af[U tn (t)] < C for t 6 [0, T e "). On the other hand, 
we have 

Af[U £ "(t)} <JV[U e "(t)-U{t)]+JV[U{t)} <C*e + C, te [0,T e "). 

The uniform bound of Af[U e ™ (t)) enables us to apply the local existence result again 
to extend the solution beyond T e ". This contradicts the definition of T e "(Ci). Thus 
the proof is complete. □ 

Remark 3.2. Corollary 13. II is an extension of the corresponding theorem in |21j . 

4. Rotational dynamics of semi-classical superfluids 

The study of superfluid dynamics in response to a rotational forcing has been the 
subject of vast experimental and theoretical work in the past years. The expectation 
value of the angular momentum, i.e. 

(4.1) m £ (t):=k ^(t, x) x 1 - ■ V x ip e (t, x) dx, 

has been mainly used to describe the dynamics [4]. In particular, a nonzero value 
of m e (t) signifies the vortex nucleation in BEC experiments. For a condensate wave 
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function of the form ip e (t : x) = a e (t, x)e^ ^ t,x ' / e we obtain 

m e (t) = - / |a e | 2 x^ • V x $ e (t, x)dx + ie a'x 1 - ■ V x a e dx 

JR d JR d 

= - I \a\ 2 x x ■ V x $(t, x) dx + 0(e). 

JR d 

In the following we denote the leading order angular momentum by 
(4.2) m(i) := - / p{t, x) x 1 - ■ v(t, x) dx , 



where, as before, (p, v) = (|a| 2 , V^). We shall also use the notation 

(d) pit) ■= \ g{x)p{t,x)dx, 

jR d 

for any smooth function g{x). 

Corollary 4.1. Let f(z) = z and impose the same assumptions as before. Then, as 
e — > it holds 

(4.3) m e (i) = m(0) + ^ ((|x| 2 ) p(t) - (\x\ 2 ) Pin ) + O(e). 
Moreover we have 

d S 

(4.4) dt 1 ^^ = ^ ' V) " (f) + 2u{^ lX2)p(t) + ° (e) ' 
2 2 

where S = ^Ji^Jl denotes the trap deformation and u)\ = \{uo\ + wf) i/ie radial 
frequency. 

In the semi-classical regime, the angular momentum expectation value is dominated 
by the classical rotational effect, in contrast to the results of [HE] (where the GPE 
is considered unsealed). The result also shows that a non- isotropic frequency for the 
trapping potential contributes to the change of the angular momentum. 

Proof. Let 

n(t) := / p(t, x)x ■ v(t, x) dx 

JR d 

and we also introduce 

X*(t):= [ \x\ 2 \r(t,x)\ 2 dx = {\x\ 2 ) p{t)+ 0(e), 

jR d 

A straightforward calculation calculation then yields 

d f 

(4.5) —m(t)=fln(t)+ p{x^-V x )Vdx. 

di JRd 

When combined with a quadratic potential of the form V = \\u ■ x\ 2 , this leads to 
the relation (|4.4[) . Also, one easily verifies that 

±(\ x \ 2 ) p(t) =2n(t). 

which, upon inserting into (|4.5[) . yields (|4.3[) after an integration w.r.t. time. □ 

We also remark that the change of n(t) also depends on m{t) as well as the total 
energy. This can be seen as follows: First we calculate 

i »(*)=-/ p(x-V x )(V(x) + f(p))dx+f pv 2 dx 

UI .lod lad 



ROTATING SUPERFLUIDS 



13 



and an integration by parts gives 

— n(t)=[ (p 2 +pv 2 )dx-[ p{x-V x )V{x)dx. 

Clt jRd J R d 

Note that for f(z) = z, the total energy (|1.6|) expressed in terms of a and v reads 

For quadratic potentials V — ^\u> ■ x\ 2 , it holds (x ■ V ' x )V(x) = 2V(x), which conse- 
quently implies 

-J- n{t) = 2 (e — Qm(t) - [ pV(x) dx] - [ p{x ■ V x )V{x) dx 
di V JR d ) JR d 

= 2{E - Qm(t)) -4 I pV(x) dx. 

JR d 

Combining these calculations yields the following closed system 

— n(t)=2(E(0)-Qm(t))-2 / \uj-x\ 2 pdx, 
dt J R d 

^-m{t)=VLn{t) + {uj 2 Xl ~Lj 2 X2 ) [ x lX2 pdx. 
< at Jm d 

Finally, we note that for an isotropic harmonic confinement, i.e. V = ^u> 2 \x\ 2 , we 
have 

m{t) - Qn(t) = 0, ^ n{t) + 2Qm(t) = 2E(0) - 2uj 2 {\x\ 2 ) p{t) . 

Since in this case, it also holds 

d _ _ „ , . 2d . , 

dt x = 2n{t) = ndt m{t) 



and we simply obtain 



m (t) = ^X(t)+m(0)-^(\x\ 2 ) piii . 



2 w " ' 2 
The obtained relations yield a closed equation for m(t), i.e. 

d 2 

^2 m (t) + ( 2fi2 + 4cj 2 )m(f) = 2QE(0) + 4?ti(0V 2 - 2Quj 2 (\x\ 2 ) Pm . 
The general solution of it can be written as 

m = d cos \/4w 2 + 2fl 2 t + C 2 sin \Aw 2 + 2fl 2 t + — ^ ~4 , 

\l 2 + 2lu 2 

where C\ and C 2 are determined by m(0) and ^m(i) = iln(0). The above 
given calculations could be used to compare numerical simulations of the full NLS 
dynamics in the spirit of [H [6] . 
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